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This  work  focuses  on  the  calculation  of  transfers  inside  a  cermet  for  solid  oxide  membrane  fuel  cells  and 
electrolyzers.  A  differential  system  of  equations  presented  in  a  previous  work  is  linearized  for  low  inlet 
current  densities  using  assumptions  that  can  be  checked  quantitatively.  By  integrating  the  linearized 
equations,  we  obtain  explicit  functions  that  allow  direct  calculation  of  the  physical  quantities  describing 
the  transfers  of  the  process  inside  the  cermet.  The  functions  show  good  agreement  with  the  values 
obtained  with  the  non-linearized  system.  In  addition,  the  model  does  not  require  any  numerical  simu¬ 
lation  to  be  solved  and  can  be  implemented  in  common  spread  sheets  fairly  accurately.  A  remarkable 
dimensionless  number,  named  A,  appears  in  the  demonstration  and  is  used  for  the  calculation  of  the 
reaction  layer  thickness  of  the  cermet,  where  99.9%  of  the  charge  transfer  occurs.  This  thickness  does  not 
depend  on  inlet  current  density  or  on  the  thickness  of  the  cermet. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Different  technologies  are  currently  being  studied  for  hydrogen 
production  with  low  carbon  dioxide  emissions.  Among  these 
technologies,  electrolysis  of  water  at  low  temperature  using  poly¬ 
mer  membranes  and  noble  metals  has  already  been  extensively 
investigated  on  an  industrial  scale.  This  electrically  driven  process 
can  also  been  carried  out  at  high  temperature  with  inorganic 
membranes.  The  latter  process  presents  several  advantages 
compared  to  low  temperature  electrolysis,  i.e.  decreased  electrical 
power  consumption,  the  electrodes  can  be  made  without  using 
expensive  noble  metals  and  the  heat  that  is  required  for  the  process 
can  be  supplied  by  nuclear  reactors,  which  limits  greenhouse  gas 
production  [1]. 


*  Corresponding  author.  Institut  Europeen  des  Membranes,  Universite  de  Mont¬ 
pellier  2,  CC  047,  2  Place  Eugene  Bataillon,  34095  Montpellier  Cedex  5,  France. 
Tel.:  +33  46714929;  fax:  +33  467149235. 

E-mail  address:  mikael.dumortier@gmail.com  (M.  Dumortier). 

0378-7753 /$  -  see  front  matter  ©  2013  Elsevier  B.V.  All  rights  reserved. 
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Among  the  high  temperature  electrolysis  processes,  solid  oxide 
technology  has  been  widely  described  in  the  literature.  The  process 
can  either  be  run  at  high  (800-1000  °C),  intermediate  (600- 
800  °C)  or  low  temperatures  (400°C-600  °C).  At  high  and  inter¬ 
mediate  temperature,  the  solid  oxide  conducts  oxygen  electrons 
and  the  process  is  less  electricity-consuming  compared  to  the  low 
temperature  process,  but  the  hydrogen  produced  needs  to  be 
separated  from  steam.  At  low  temperature,  the  solid  oxide  can 
conduct  protons  2],  which  means  the  hydrogen  produced  at  the 
cathode  is  pure  and  a  subsequent  separation  step  can  be  avoided 
(Fig.  1).  In  addition,  the  lower  temperature  process  contributes  to 
the  long-term  resistance  of  materials. 

In  high  temperature  electrolysis,  the  cell  is  an  electrochemical 
membrane  reactor  composed  of  two  channels  separated  by  a  dense 
solid  oxide  layer  between  two  porous  electrodes.  In  proton  con¬ 
ducting  cells,  steam  is  fed  at  the  anode  and  oxidized  in  the  cermet, 
which  is  a  mixed  media  made  of  metal  and  ceramic.  Mathematical 
modelling  of  such  devices  helps  understanding  the  phenomena 
inside  the  cells  and  makes  it  easier  to  predict  their  operation. 


704 


M.  Dumortier  et  al.  /  Journal  of  Power  Sources  248  (2014)  703-713 


Nomenclature 

molar  fraction  for  species  i  (dimensionless) 

A 

transfer  current  density  distribution  number 

^ox>  ^red 

molar  fraction  for  the  oxidizer  and  the  reducer 
(dimensionless) 

(dimensionless) 

a 

charge  transfer  coefficient  (dimensionless) 

Cj 

total  concentration  (mol  m-3) 

£i 

volumetric  fraction  of  the  phase  i  (dimensionless) 

Di 

diffusion  coefficient  for  species  i  (m2  s-1) 

V 

overpotential  (V) 

Deff 

ui,  Kn 

effective  Knudsen  diffusion  coefficient  for  species  i 

Vi 

stoichiometric  coefficient  for  species  i  (dimensionless) 

(m2  s-1) 

Pg 

gas  density  (kg  m-3) 

Deff 

y 

effective  binary  diffusion  coefficient  for  species  i  and  j 

Vi 

electrical  conductivity  for  phase  i  (S  m-1) 

(m2  s-1) 

(pi 

charge  potential  inside  phase  i  (V) 

~ex 

unitary  vector  directed  to  the  thickness  of  the 

Vred 

ceramic  conductivity  fraction  (dimensionless) 

F 

electrode  (dimensionless) 

Faraday  constant  (C  mol-1) 

Numerical  analysis 

I 

inlet  intensity  (A) 

V 

del  operator 

ir 

volumetric  transfer  current  (A  m-3) 

A 

change  of  variable  quantity 

i° 

tr) 

exchange  current  density  (A  m-3) 

d 

differential  operator 

h 

charge  flux  density  inside  phase  i  (A  m-2) 

0 

partial  differential  operator 

Jo 

k 

inlet  electrical  current  density  (A  m-2) 

charge  of  the  ion  conducted  by  the  ceramic  phase 

Indexes  and  subscripts 

I< 

electrochemical  transfer  number  (dimensionless) 

e 

variable  attached  to  the  metallic  phase  or  to  charge 

L 

electrode  thickness  (m) 

carried  by  electrons 

f'reac 

reaction  layer  thickness  (m) 

ext 

used  for  surfaces  included  in  the  boundary  of  the 

Mi 

molar  weight  for  species  i  (kg  mol-1) 

control  volume 

Mi 

molar  weight  of  the  transported  ion  (kg  mol-1) 

eff 

effective  value  of  a  variable 

nT 

molar  flux  density  for  species  i  (mol  s-1  m-2) 

g 

variable  attached  to  the  gas  phase 

P 

total  pressure  (Pa) 

io 

variable  attached  to  the  ceramic  phase  or  to  charge 

Pe 

Peclet  number  (dimensionless) 

carried  by  ions 

Q 

volumetric  flow  rate  (m3  s-1) 

int 

used  for  surfaces  inside  the  control  volume 

R 

ideal  gas  constant  (J  mol-1  K-1) 

TP 

variable  attached  to  the  triple  phase  boundary 

n 

grain  radius  inside  phase  i  (m) 

* 

dimensionless  quantity 

rr 

total  volumetric  molar  transfer  rate  (mol  s-1  m-3) 

0 

reference  quantity 

Si 

T 

surface  fraction  for  phase  i  (dimensionless) 
local  uniform  temperature  (I<) 

Abbreviations 

u 

velocity  vector  (ms-1) 

PDE 

partial  differential  equation 

X 

arbitrary  physical  quantity 

SOEC 

solid  oxide  electrolysis  cell 

x+ 

reference  added  concentration  factor  (dimensionless) 

SOFC 

solid  oxide  fuel  cell 

X 

position  across  the  electrode  (m) 

Fortunately  for  modellers,  the  equations  describing  transport 
phenomena  inside  solid  oxide  fuel  cells  (SOFC)  and  solid  oxide 
electrolysis  cell  (SOEC)  are  the  same,  as  one  technology  is  the 
reverse  functioning  of  the  other.  Therefore  many  numerical  models 
are  valid  for  both  technologies.  While  complex  3D  models  have 
been  developed  to  describe  transfers  in  high  temperature  fuel  cells 
or  electrolyzers  [3-15],  only  few  have  integrated  models  available 
for  the  non-code  user.  In  2007,  Meng  Ni  proposed  a  ID  analytical 
model  to  predict  the  potential  inside  a  solid  oxide  electrolysis  cell 
(SOEC)  without  numerical  solving  [16  .  The  model  is  based  on  the 
assumption  that  the  electrochemical  reaction  is  confined  to  a  thin 
layer;  in  a  previous  work,  our  group  demonstrated  that  the 
assumption  is  only  valid  under  certain  conditions  [17],  i.e.  low  ionic 
conductivity  of  ceramics,  no  limitation  in  the  diffusion  of  species 
inside  the  electrode,  metal  grain  radius  not  too  large  compared  to 
the  ceramic  grain  radius.  Flowever,  these  statements  are  only 
qualitative.  The  objective  of  this  paper  is  to  provide  a  model 
featuring  analytical  expressions  of  the  physical  quantities  inside  a 
cermet  electrode  used  in  solid  oxide  cells,  assuming  that  the  elec¬ 
trochemical  reaction  can  occur  anywhere  in  the  cermet.  The  model 
is  derived  from  a  comprehensive  3D  model  which  has  been 
simplified  with  two  major  hypotheses:  unidirectional  transport 
phenomena  and  linearized  Butler— Volmer  equation.  The  frames  in 
which  those  assumptions  are  valid  are  also  developed  in  this  study. 


The  model  does  not  require  any  numerical  simulation  to  be  solved 
and  can  be  implemented  in  common  spread  sheets  fairly  accu¬ 
rately.  In  order  to  be  the  most  general  possible,  the  model  equations 
are  written  in  dimensionless  form.  The  undimensionalization  of  the 
governing  equations  for  transfer  inside  a  fuel  cell  has  already  been 
investigated  by  Gyenge  in  2005  18].  Flowever  to  our  knowledge,  no 
integration  of  a  linearized  system  of  equations  has  ever  been 
investigated  to  predict  the  value  of  the  physical  quantities  along  a 
porous  electrode. 

2.  Modelling 

In  this  study,  the  electrodes  considered  are  cermets,  which  are 
among  the  most  interesting  candidates  as  electrode  material  for 
high-temperature  steam  electrolysis  [2].  A  cermet  is  a  mixed 
porous  media  comprising  three  phases:  a  metallic  phase  which 
conducts  electrons,  a  ceramic  phase  which  conducts  protons  and  a 
gas  phase  constituted  by  the  porosity  where  reactants  and  products 
of  the  electrochemical  reaction  are  conveyed  (Fig.  2). 

While  charge  is  transported  by  migration  in  solid  phases, 
chemical  species  are  transported  in  gas  by  diffusion  and  convec¬ 
tion.  These  transport  phenomena  inside  electrodes  have  been 
mathematically  modelled  in  a  previous  work  17]  using  continuity 
equations  obtained  by  applying  the  volume  averaging  method  [19] 
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Fig.  1.  Principle  of  proton  and  anion  conducting  SOECs. 


to  a  representative  volume  of  electrode.  As  explained  earlier,  this 
paper  focuses  on  the  establishment  of  a  ID  model  derived  from 
these  continuity  equations. 


2  A.  Assumption  1:  one  dimensional  model 

The  transport  phenomena  are  represented  by  partial  differential 
equations  (PDE).  A  PDE  for  the  function  u(x,y,z)  is  an  equation  of  the 
form  : 


d2U 


9  2u 


dU  dU 

x’^’z,u,dx’dy"  'dydz' ' '  dz2 


=  0 


(1) 


X  =  Xj  X  =  x2 

X*  =  0  X*=l 


t 

Electric  current 

Gas  flow 

Channel 

Electrode 

Electrolyte 

> 

y 

rrx 

\  x 

_ - 

X 


Fig.  3.  Formal  representation  of  the  electrolyser. 


be  perpendicular  to  the  electrode.  As  a  result,  the  divergence 
operator  and  the  gradient  operator  can  be  simplified  as  follows: 


(3) 


Different  conditions  are  required  to  validate  this  assumption. 
First  of  all,  the  electric  current  must  be  quite  homogeneous  along 
the  y  direction  of  the  cermet.  This  statement  presumes  that  the 
distribution  of  the  metallic  and  ceramic  phases  along  the  y  direc¬ 
tion  is  homogeneous.  The  quantities  can  however  vary  across  the  x 
direction  of  the  electrode. 

Secondly,  during  the  electrochemical  process,  the  mixture  in  the 
anode  channel  is  enriched  in  oxygen  or  hydrogen  depending  on  the 
type  of  membrane  reactor,  i.e.  a  fuel  cell  or  electrolyser.  If  the  flow  is 
crosswise  to  the  electrode,  the  variation  of  molar  fractions  in  the 
channel  in  the  y  direction  can  have  an  impact  on  the  transport 
phenomena  in  the  cermet.  With  a  simple  analysis,  the  changes  in 
molar  fractions  inside  the  channel  can  be  established  through 
Equation  (4).  In  this  work  we  believe  that  the  1 D  assumption  can  be 
reasonably  considered  valid  if  the  variations  are  lower  than  10%. 


A  PDE  often  requires  numerical  simulation  in  order  to  be  solved. 
However,  using  this  effective  calculation  strategy  gives  the  values  of 
u(x,y,z)  but  not  the  expression  of  the  function  itself.  One  aim  of  this 
work  was  to  find  assumptions  that  can  be  realistically  considered 
regarding  the  behaviour  of  the  system  to  reduce  those  PDE  to  a 
simpler  form. 


.  RTI  Vi  n . 


(4) 


Taking  into  consideration  this  assumption,  the  analytical  ex¬ 
pressions  of  the  variables  of  the  system  can  be  obtained. 


F(x,u)  =  0  (2) 

To  obtain  the  functions,  we  assumed  that  the  different  quanti¬ 
ties  vary  along  a  single  direction  x,  oriented  by  the  electric  current 
path  across  the  electrode  (Fig  3).  Here,  the  direction  is  assumed  to 


2.2.  Current  densities  and  potentials 

In  a  cermet,  electrons  are  transported  by  the  metal  and  ions  by 
the  ceramic.  The  behaviour  of  this  transport  is  represented  by 
Ohm’s  law: 


Ceramic  phase 


Metallic  phase 


Gas  phase 


Triple  phase  boundary 


Fig.  2.  The  cermet  medium  and  its  different  phases. 


A 

For  the  ceramic  :  jio  =  -<7io-^  (5) 

For  the  metal :  je  =  (6) 

9x 

The  transport  number  of  charged  species  inside  the  ceramic  is 
supposed  equal  to  1,  meaning  that  only  one  ion  is  conveyed  in  the 
ceramic,  hence  mixed-conducting  ceramics  are  not  studied  here. 
During  the  electrochemical  reaction,  charge  is  transferred  from  the 
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metal  to  the  ceramic.  The  rate  of  this  transfer  is  calculated  using  the 
Butler— Volmer  equation  [17]: 


ir  = 


P^exp 

LX?ed 


(7) 


helps  in  the  comparison  of  the  different  terms,  and  can  also  be  used 
to  recover  characteristic  properties  of  the  system  such  as  length, 
time  constants  and  dimensionless  numbers.  The  new  variables  are 
expressed  as  the  ratio  of  the  variables  to  a  relevant  order  of 
magnitude,  as  in  this  example: 


The  overpotential,  p,  corresponds  to  the  difference  of  electrical 
state  between  the  two  solid  phases: 

V  =  (<Pe  -  <Pio )  -  (<Pe  -  K)  (8) 

The  reference  electric  potential,  0°,  is  calculated  with  the  Gibbs 
energy  of  the  electrochemical  reaction  and  the  reference  ionic 
potential,  0^,  is  set  to  0. 

As  charge  passes  from  a  medium  to  another,  the  value  of  the 
currents  varies  with  the  intensity  of  the  charge  transfer  density,  ir. 
The  variables  are  interconnected  through  the  charge  conservation 
expressions: 


X*  =  X/X0  (16) 

The  reference  value  depends  on  the  parameters  of  the  process, 
and  is  not  affected  by  derivation: 


r*  _  ax^  _  xq  ox 
“  ax*  ”  x0  ax 


(17) 


Reference  values  are  obvious  for  position  and  current  densities: 


x/l 

(18) 

ir /i? 

(19) 

For  the  ceramic  :  ^O'ioSio)  =  I’ehtw  (9) 

For  the  metal  :  f(jese)  =  -veirejp  (10) 

OX 

Finally,  current  conservation  inside  the  phases  is  expressed  by 
the  equation  below: 

7io5io  +  7e5e  —  Jo  ( 1 1 ) 

2.3.  Assumption  2:  low  currents  and  overpotentials 


Je  =  JeSe/Jo  (20) 

j[o  —  Jio5io//o  (21) 

The  reference  values  for  overpotentials  are  deduced  from  But- 
ler-Volmer  equation: 

V*  =  V/Vo  (22) 

Vo  =  RT/aveF  (23) 


The  dimensionless  forms  of  Equations  (5) — (11 )  are  written  as 
follows: 


This  equation  system  cannot  be  solved  analytically  because  of 
the  non-linearity  of  the  Butler-Volmer  equation.  However,  the 
exponential  terms  of  the  equation  can  be  linearized  by  a  first  order 
Taylor  development  at  0  as  follows: 


)\o  —  (Siof?io/<7o)dx*4>io 

(24) 

Je  =  -(Se<T  e/<7o)dx*0e 

(25) 

ex  =  1  +  x  +  o(x)  (12) 

This  assumption  is  valid  for  low  overpotential  or  low  current 


densities: 

aveF 

~Kfr]  <  1 

(13) 

Jo  -j 

l?eJPL  « 

(14) 

Note  that  Equation  (13)  states  that  a  current  is  considered  to  be 
low  when  it  is  much  smaller  than  the  exchange  current  density  of 
the  electrode,  i.e.  when  the  system  is  close  to  the  thermodynamical 
equilibrium.  Using  Equation  (4),  the  variation  of  the  molar  fractions 
of  species  with  x  is  also  supposed  to  be  negligible  compared  to  the 
variation  of  the  overpotential.  Therefore  Equation  (7)  reduces  to: 


fr  =  xV 

(26) 

dWe  =  -«r 

(27) 

dx*Jio  =  Kir 

(28) 

.*  .*  . 

Je  +7  io  1 

(29) 

with 

X+  =  ^red/x?ed +xoxAox 

(30) 

=  JqL/ Vo 

(31) 

ir  = 


0aVeF 

r  RT 


(15) 


2.4.  Undimensionalization  of  the  system 

Once  these  assumptions  are  made,  the  equations  can  be  written 
in  dimensionless  form.  The  process  of  undimensionalization  con¬ 
sists  in  making  a  change  of  variables  in  equations  so  that  the  new 
variables  and  their  derivatives  are  of  the  same  order  of  magnitude. 
It  facilitates  the  manipulation  of  equations  during  calculations, 


I<  =  Lvetfejp/Jo  (32) 

Coupling  the  Butler-Volmer  equation  (26)  with  the  continuity 
equation  for  charge  ((27)-(28))  and  current  conservation  (29)  gives 
a  second  order  differential  equation  verified  by  the  ionic  charge 


current: 

till- Ml  =  *»* 

(33) 

with  _  si0<ri0 

Fed  Sio^io  +  5e^e 

(34) 
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A  ,■  (np+%)  (- !— +-L)  (35) 

RT  Wed  ^ox/  V’io'Ao 

The  boundary  conditions  for  the  ionic  current  are  given  by: 

No  ionic  current  at  the  electrode/channel  interface  :j*Q( 0)  =  0 

(36) 

At  the  electrolyte/electrode  interface  the  current  is  purely  ionic 
=  1 

(37) 

The  solution  of  this  differential  system  composed  of  Equations 

(33),  (36)  and  (37)  is: 


=  ffred(l  +/+  (x*)  -/  (x*)) 


with 


f+(x *)  =  De  exp(v^x*) 
/_(x*)  =  (De  +  l)exp(-v(Ax*) 


De  = 


|^  +  exp(-vA 
exp^v^)  -exp(— v^) 


(38) 

(39) 

(40) 

(41) 


From  Equations  (29)  and  (38),  the  expression  of  dimensionless 
electrical  current  is  obtained: 

fe(x)  =  1-^(1  +/+(*’)  -/-(*’))  (42) 

By  derivation  of  this  quantity,  and  coupling  it  with  Equation 
(27),  the  expression  of  the  exchange  current  density  in  dimen¬ 
sionless  form  is  given  by: 


lr  = 


Va 

K 


<Aed  (/+ (**)+/  (**)) 


(43) 


Elence,  the  dimensionless  overpotential  using  Equation  (26)  can 
be  written  as  follows: 


v  =^ared(r(x)  +/-(x*)) 


2.5.  Mass  transfer 


(44) 


The  chemical  species  needed  for  the  electrochemical  reactions 
are  transported  to  the  triple  phase  boundary  in  the  porous  elec¬ 
trode.  The  model  chosen  for  mass  transport  inside  solid  oxide  cells 
is  the  Dusty  Gas  Model,  as  recommended  by  Suwanwarangkul  et  al. 
[4].  For  two  species,  we  have  shown  in  our  previous  work  [17]  that 
the  Dusty  Gas  model  can  be  expressed  as  follows: 


N,  =  -D~Ct  +  XjU 


with 


D,  = 


DeU  Deff 

i,I<n  ij 


(45) 


(46) 


Df+xJDl  L*XiD!l 

Reducing  this  equation  to  its  dimensionless  form  gives: 

N-  =  -dx*Xi/Pei  +  u*Xj  (47) 


With 


N  = 


Nix$  g 

N? 

(48) 

XiCTu/N° 

(49) 

kn 

F  ve 

(50) 

LNf  =  RTJ0L1 >± 

(51) 

CjDj  P  FDj  ve 

pei  = 


Pei  is  the  Peclet  number  for  species  i.  It  compares  advection  and 
diffusion  effects  in  mass  transfer.  If  Pez  <<  1,  then  transport  of 
chemical  species  inside  the  medium  occurs  by  diffusion  while  if 
Pei  »  1,  then  mass  transfer  is  governed  by  advection. 

In  this  work,  the  Peclet  number  is  supposed  to  be  very  small  in 
the  electrodes.  In  the  problem  presented  in  Section  3,  the  Peclet 
number  of  oxygen  inside  the  porous  anode  for  a  proton  con¬ 
ducting  SOEC  is  equal  to  Pe 02  =  2.8  x  10-3  (Pe h2  =  1.7  x  10~2  in 
the  cathode).  This  reasonably  implies  that  the  velocity  inside  the 
electrode  is  low  and  that  the  chemical  species  are  transported  by 
diffusion  inside  the  cermet.  Elence,  Equation  (47)  can  be  written  as 
follows: 


n;  =  -dx*Xi/pei 


(52) 


In  addition,  chemical  species  are  either  consumed  or  produced 
inside  the  porous  electrodes.  Therefore  their  fluxes  depend  on  the 
intensity  of  the  electrochemical  reaction.  The  interconnection  be¬ 
tween  the  phenomena  is  given  by  mass  conservation  for  the 
species: 


fa(SSNi)  =  vii  rFeTP 


(53) 


the  dimensionless  form  of  the  equation  above  can  be  written  as: 

dx*N;  =  Ki*  (54) 

Similarities  can  be  noted  between  dimensionless  Equations  (28) 
and  (54).  This  shows  that  charge  and  mass  conservation  exhibit  the 
same  physical  behaviour.  This  is  consistent  with  the  fact  that  the 
transport  of  those  quantities  are  driven  and  bounded  together  by 
the  electrochemical  transfer. 

The  molar  flux  density,  Nz,  is  obtained  through  integration 
of  (54).  Hence,  boundary  conditions  are  required.  At  the  elec¬ 
trode/electrolyte  interface,  the  molar  flux  of  species  is  equal  to 
zero  as  they  cannot  cross  the  electrolyte.  This  gives  the 
condition: 


Mj  (1)  =  0 


(55) 


At  the  electrode/channel  interface,  the  molar  flux  of  the  species  i 
is  equal  to  the  total  source  of  the  species  in  the  cermet. 

The  mass  transport  variables  can  be  then  calculated  from  the 
integration  of  Equations  (52)  and  (54): 


N,*(x)  =  -fe(x) 

Axf  (x*)  =  Pei  (x*  -  aTed  (x*  +  (/+  (x*) 
+/-(**)  -2De-l)/vft)) 


(56) 

(57) 
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2.6.  Additional  assumption  and  physical  meaning  of  the  A  number 

In  most  cases,  electrical  conductivity  of  metals  ranges  between 
103  et  105  S  m_1.  These  values  are  at  least  100  times  higher  than  the 
ionic  conductivity  of  ceramics  [20—22].  Hence  seae  »  Siocrio.  When 
this  assumption  is  injected  into  the  previous  equations,  they 
simplify  into  much  simpler  expressions  and  exhibit  a  strong 
dependence  on  the  A  number: 


Fig.  4.  Dimensionless  transfer  current  density  for  several  values  of  A.  X-axis:  dimen¬ 
sionless  position  in  the  cermet.  Y-axis:  dimensionless  transfer  current  density. 


Ax, 


1  ch| 

[VAx*) 

Va 

sh  | 

)  > 

(61) 


The  equations  show  that  the  A  number  is  a  key  dimensionless 
number  describing  the  distribution  of  the  electrochemical  reaction 
inside  the  electrode.  Equation  (35)  shows  which  process  and 
structural  parameters  are  involved  into  the  calculation  of  A , 
therefore  parametric  analysis  to  study  the  influence  of  parameters 
on  the  quantities  related  to  the  electrochemical  reaction  can  easily 
be  made.  To  our  knowledge,  this  dimensionless  number  has  never 
been  described  in  the  literature.  The  number  compares  two 
limiting  effects  that  determine  the  location  of  charge  transfer: 


•  The  cermet  resistance  to  charge  transfer,  which  depends  on  the 
microstructure,  gas  composition,  electrode  thickness  and  ex¬ 
change  current  density.  Electrochemical  transfer  conductivity  of 
the  electrode,  <7th,  is  the  characteristic  conductivity  of  the  elec¬ 
trode  for  the  electrochemical  transfer. 

•  Charge  conduction  resistance  inside  the  metal  and  the  ceramic. 
When  charge  is  transported  inside  the  conductors,  it  generates 
energy  loss  because  of  resistance  effects.  The  effects  are  repre¬ 
sented  in  the  A  number  by  the  conductivity  oohm.  which  is  the 
harmonic  mean  of  the  conductivities  of  the  two  charge 
conductors. 


A  =  ffth 

(62) 

^OHM 

L2v2ei°eTPaF  (xTeA  :  Xox\ 
th  RT  Wed  X°oJ 

(63) 

1  11 

(64) 

0OHM  Sio^io  5e^e 

As  shown  in  Fig  4,  when  charge  conduction  transport  is  limiting, 
then  (7th  >>  0ohm  and  A  »  1,  and  the  reaction  occurs  mostly  at  the 
anode/electrolyte  interface.  This  is  consistent  with  our  previous 
observations  17],  in  which  we  noted  that  when  ceramic  conduc¬ 
tivity  is  low,  charge  transfer  occurs  close  to  the  electrolyte  in  order 
to  minimize  the  path  of  protons  through  the  ceramic. 


When  electrochemical  transfer  is  limiting,  i.e  crt h  <<  oohm  and 
A  «  1,  the  reaction  tends  to  follow  the  distribution  of  at h  across  the 
electrode.  Here,  ath  is  homogenous  across  the  electrode,  but  the 
profile  would  be  different,  for  example,  in  the  case  of  graded 
composition  electrodes  [6,23-25]. 

The  distributions  of  transport  fluxes  and  transfer  fluxes  across 
the  electrodes  only  depend  on  A.  Equations  (62)-(64)  show  how 
the  parameters  of  the  process  influence  those  distributions.  For 
example,  inlet  current  density  does  not  have  any  effect  on  the 
distributions,  contrary  to  the  number  of  electrons  exchanged  in  the 
electrochemical  transfer.  From  the  expression  of  A,  we  can  conclude 
that  distribution  of  the  electrochemical  reaction  across  the  elec¬ 
trode  depends  mainly  on  the  reaction  involved,  on  conductivities 
and  on  the  electrode  microstructure,  which  again,  is  consistent 
with  previous  observations  [17]. 

2.7.  Dimensional  expressions 

Dimensional  expressions  are  required  to  find  real  values  on  a 
real  electrode.  Dimensional  expressions  are  obtained  from  corre¬ 
sponding  dimensionless  expressions  with  the  relationships  pre¬ 
sented  below: 

x  =  X*(X!  -x2)  +Xj  (65) 

X(x)  =  X0X*  (x*)  =  X0X*((x  -  x,)/{x,  -  x2))  (66) 

Xi  and  X2  are  the  x  coordinates  of  the  channel/electrode  interface 
and  electrode/electrolyte  interface  respectively.  The  expressions 
are  used  for  the  comparison  between  the  ID  model  presented  in 
Section  2.3  and  the  3D  complete  model  presented  in  Section  2.1. 

3.  Comparison  between  ID  and  3D  models 

3.1.  Calculation  set  up 

In  a  previous  study  [17],  the  3D  model  results  were  presented  as 
a  planar  co-flow  proton  conducting  SOEC  (see  Fig.  3)  using  the 
COMSOL  Multyphysics  4.2  software.  The  study  showed  that  dis¬ 
tribution  of  the  electrochemical  reaction  across  the  cermet  depends 
on  several  parameters  such  as  ionic  conductivity  and  electrode 
microstructure.  These  results  will  be  used  as  a  basis  to  test  the 
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simplified  model  presented  here,  considering  the  same  reactor 
dimensions  (Table  1)  and  boundary  conditions  (Table  2),  parame¬ 
ters  and  reference  values  used  in  our  previous  work.  Unless  another 
value  for  the  parameters  is  specified,  the  default  values  taken  for 
any  calculation  are  those  presented  in  Table  3.  In  order  to  facilitate 
reading,  the  model  developed  in  the  previous  work  will  be  referred 
to  as  the  “3D  model”  whereas  the  model  developed  in  this  work 
will  be  referred  to  as  the  “ID  model”.  The  formulae  presented  in 
Section  2.3  are  expressed  in  dimensional  form  according  to  the 
geometry  presented  in  Fig.  5. 

Transfer  current  density,  iV^tbp,  is  proportional  to  the  charge 
transferred  at  the  triple  phase  boundary  and  to  the  production  and 
consumption  rate  of  the  chemical  species.  Therefore,  it  is  connected 
to  the  location  where  the  electrochemical  reaction  occurs  in  the 
electrode.  The  integral  of  this  quantity  is  fixed  by  the  inlet  current 
density  value.  Therefore  other  parameters  can  only  have  an  effect 
on  its  distribution  across  the  electrode.  Transfer  current  density  is 
thus  the  best  suited  quantity  for  this  study  and  for  a  parametric 
analysis.  The  two  models  are  therefore  compared  in  terms  of  their 
prediction  of  transfer  current  density,  and  also  in  terms  of  their 
prediction  of  the  overpotential,  77,  and  the  oxygen  molar  fraction, 
X02.  With  the  latter  two  quantities,  the  differences  between  both 
models  should  appear  more  evidently,  as  they  are  less  constrained 
than  i>TP. 

The  specific  reactive  volume  etp  can  either  be  determined 
experimentally  of  calculated.  Several  models  can  be  found  in 
literature  to  relate  this  volume  to  the  structural  parameters  of  the 
electrode  [26,27].  Correlations  for  the  calculation  of  binary  and 
Knudsen  diffusion  coefficient  can  be  obtained  from  the  well-known 
works  of  Bird  et  al.  and  Poling  et  al.  [28,29].  Before  running  the 
calculations  for  the  comparison,  one  has  to  check  whether  the  as¬ 
sumptions  for  the  1 D  model  are  satisfied.  As  the  maximum  current 
density  in  this  study  is  quite  high  (10,000  A  rrr2),  the  linearization 
of  the  Butler  Volmer  equation  in  the  ID  model  could  be  questioned. 
Yet,  from  the  values  for  the  exchange  current  density  and  the 
specific  reactive  volume  from  the  previous  study,  the  low  current 
assumption  presented  in  Equation  (13)  is  verified  for  the  reference 
values: 


Jo 

i^ejpL 


0.15 


(67) 


The  length  of  the  cell  (1  mm)  is  below  the  usual  length  found  in 
literature.  In  fact,  the  priority  of  this  comparison  is  not  to  represent 
real  system  but  to  test  the  ID  model  in  a  cell  where  the  assump¬ 
tions  explained  above  are  satisfied.  Therefore,  the  length  of  the  cell 
was  purposely  kept  low  to  reduce  calculation  cost  and  satisfy  the 
ID  assumption  presented  in  Equation  (12).  From  the  parameters  of 
Tables  1  and  2,  the  ID  assumption  is  verified  for  every  species: 


Table  2 

Boundary  conditions. 


Boundary 

Species  transport 

Momentum 

transport 

Charge  transport 

1 

Zero  flux 

Wall  —  no  slip 

- 

2 

Inflow  —  concentration 

Velocity  inlet 

- 

3 

Outflow 

Outlet  -  pressure 

- 

4 

Zero  flux 

Wall  —  no  slip 

- 

5 

- 

- 

Ionic:  zero  flux 

6,8,10,11, 

Zero  flux 

Wall  —  no  slip 

electronic:  flux 

Zero  flux 

13,14 

7,9,12,15 

Zero  flux 

Wall  —  no  slip 

16 

Outflow 

Outlet  -  pressure 

- 

17 

- 

- 

Ionic:zero  flux 

18 

Zero  flux 

Wall  —  no  slip 

electronic:  Flux 

19 

Inflow  —  concentration 

Velocity  inlet 

- 

20 

Zero  flux 

Wall  —  no  slip 

- 

A  higher  length  results  in  higher  intensities  I  for  the  same  inlet 
current  density. 

3.2.  Parametric  results 

The  two  models  were  firstly  compared  with  the  reference  values 
given  in  Table  3  and  then  for  different  conductivity  ratios.  In  both 
cases,  comparison  between  the  ID  and  the  3D  models  showed  a 
perfect  match  for  TP  (see  Figs.  6  and  7). 

A  good  agreement  was  also  achieved  for  prediction  of  the  oxy¬ 
gen  molar  fraction,  x02,  and  the  overpotential,  r\.  The  results  pre¬ 
sented  in  Figs.  8  and  9  show  that  the  ID  model  tends  to  slightly 
underestimate  the  value  of  the  overpotential  when  its  value  is  low. 
This  corresponds  to  the  effect  of  the  term  (xred/x°ed  -  x0X  AqX) 
which  was  neglected  in  Section  2.2  in  order  to  analytically  solve 
Equation  (33).  The  influence  of  this  term  becomes  important  when 
the  overpotential  is  low.  The  values  of  oxygen  molar  fraction  across 
the  cermet  are  plotted  in  Figs.  10  and  11.  Again  the  ID  model 
perfectly  matches  the  3D  model  except  when  the  ratio  between  the 
charge  conductivities  of  the  ceramic  and  the  metal  is  low. 

In  addition  to  comparing  the  two  models,  the  results  show  that 
when  the  ratio  between  the  conductivities  is  high,  the  electro¬ 
chemical  reaction  occurs  in  an  area  close  to  the  electrolyte.  This 
phenomenon  has  been  explained  in  our  previous  work  [17]:  the 
profile  adapts  in  order  to  minimize  the  path  of  protons  through  the 
ceramic  to  lower  ohmic  losses. 

The  results  obtained  on  comparison  of  the  two  models  for  other 
parameters  (grain  ratio  radii,  initial  oxygen  molar  fraction,  inlet 
current  density)  are  not  presented  here  as  they  can  be  described  in 


RTI  Vi 

--^<0.019 

PQF  ve- 

(68) 

Table  1 

Reactor  dimensions. 

Parameter 

Value 

Anodic  channel  thickness 

4  x  KT4  m 

Anodic  cermet  thickness 

1.5  x  1CT4  m— 5  x  1CT4  m 

Electrolyte  thickness 

1  x  10-4  m 

Cathodic  cermet  thickness 

1.5  x  1(T4  m— 5  x  1CT4  m 

Cathodic  channel  thickness 

4  x  1CT4  m 

Anodic  cermet  length 

1(T3  m 

Electrolyte  length 

1(T3  m 

Cathodic  cermet  length 

1(T3  m 

Anodic  channel  length 

1.8  x  1CT3  m 

Cathodic  channel  length 

1.8  x  10"3  m 

Table  3 

Values  of  variable  parameters  and  inlet  parameters  for  the  reference  calculation. 


Parameter 

Reference 

value 

Unit 

Other  values  tested 

Inlet  current  density 

10,000 

A  ITT2 

1-1000 

Ionic  conductivity 

1 

S  m  1 

- 

(ceramic) 

Electronic  conductivity 

1000 

S  m  1 

1-100,000 

(metal) 

Particle  radius  (metal) 

10"6 

m 

5  x  10-7— 5  x  10"6 
-2.5  x  10"5 

Particle  radius  (ceramic) 

5  x  10~7 

m 

- 

Initial  oxygen  molar 

0.01 

- 

0.1,03,  0.5,  0.9 

fraction 

Operating  pressure 

106 

Pa 

- 

Operating  temperature 

873 

K 

- 
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Fig.  5.  Geometry  used  for  the  calculations  and  boundary  conditions. 


similar  terms:  perfect  match  for  i>Tp  and  x02,  small  difference  for  77. 
In  addition,  we  have  noted  that  the  difference  between  the  models 
for  77  decreases  with  lower  inlet  current  density.  However,  addi¬ 
tional  calculations  showed  that  the  ID  model  deviates  from  the  3D 
model  when  the  low  current  assumption  (Equation  (13))  is  no 
longer  valid.  Therefore  we  can  highlight  the  fact  that  the  frame  of 
the  1 D  model  is  restricted  to  low  currents. 


x  (10-4  m) 

Fig.  7.  Comparison  between  the  two  models  on  transfer  current  density  along  the 
anodic  cermet  for  different  conductivities  ratios.  X-axis:  position  in  the  anodic  cermet 
(on  the  left  the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  transfer  current 
density  multiplied  by  the  specific  reactive  volume. 

3.3.  Reaction  layer  thickness 

One  of  the  most  interesting  results  of  this  work  comes  from  the 
parametric  study  through  the  electrode  thickness  (see  Fig  12). 
Firstly,  we  noticed  that  the  distribution  of  transfer  current  density 
is  independent  of  the  thickness  of  the  electrode.  This  means  that 
the  part  of  the  electrode  that  contributes  to  the  reaction  has  a 
characteristic  length  which  depends  on  several  other  parameters  of 
the  process.  The  objective  of  the  demonstration  that  follows  is  to 
find  an  analytical  expression  for  this  characteristic  length. 

There  are  different  methods  to  obtain  the  expression  of  the  re¬ 
action  layer  thickness.  The  one  presented  here  is  the  most  rigorous. 
We  started  with  the  fact  that  the  total  charge  exchanged  during  the 
electrochemical  reaction  is  determined  by  the  value  of  the  inlet 
current  density.  This  charge  is  exchanged  from  one  conductor  to 
another  all  over  the  cermet  through  the  electrochemical  reaction. 
Hence: 

i 

Jo  =  J  ve  Sjpirdx  (69) 

0 


x  (10"4m) 

Fig.  6.  Comparison  between  the  two  models  on  transfer  current  density  along  the 
anodic  cermet  for  the  reference  case.  X-axis:  position  in  the  anodic  cermet  (on  the  left 
the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  transfer  current  density 
multiplied  by  the  specific  reactive  volume. 


x  (10-4  m) 

Fig.  8.  Comparison  between  the  two  models  on  overpotential  along  the  anodic  cermet 
for  the  reference  case.  X-axis:  position  in  the  anodic  cermet  (on  the  left  the  anodic 
channel,  on  the  right  the  electrolyte).  Y-axis:  overpotential. 
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Fig.  9.  Comparison  between  the  two  models  on  overpotential  along  the  anodic  cermet 
for  different  conductivities  ratios.  X-axis:  position  in  the  anodic  cermet  (on  the  left  the 
anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  overpotential. 

We  can  define  the  reaction  layer  thickness,  Lreao  as  the  portion 
of  electrode  in  which  99.9%  of  this  charge  is  exchanged. 

L 

999  r 

^ Jo  =  J  ^plrdx  (70) 

L—Lye  ac 

The  assumption  presented  in  Section  2.4  on  conductivities  was 
then  considered.  Moreover,  defining  a  reaction  layer  supposes  that 
the  electrochemical  reaction  is  confined  to  an  area  close  to  the 
electrolyte,  which  is  valid  for  values  of  A  »  1.  Therefore  the  hy¬ 
perbolic  functions  can  be  approximated  as  exponential  functions: 

chx-shx-0^3*  (71) 

using  Equation  (71)  for  the  integration  of  Equation  (70)  : 

1000  =  exp(_v^*reac)  (72) 

We  can  now  use  a  normal  logarithm: 

W  =  7/[  (73) 


x(10'4m) 


Fig.  10.  Comparison  between  the  two  models  on  the  molar  fraction  of  oxygen  along 
the  anodic  cermet  for  the  reference  case.  X-axis:  position  in  the  anodic  cermet  (on  the 
left  the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  molar  fraction  of  oxygen. 


Fig.  11.  Comparison  between  the  two  models  on  the  molar  fraction  of  oxygen  along 
the  anodic  cermet  for  different  conductivities  ratios.  X-axis:  position  in  the  anodic 
cermet  (on  the  left  the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  molar 
fraction  of  oxygen. 


This  expression  is  valid  for  A  >  50.  If  A  <  50,  then  Lreac  >  L  and 
defining  a  reaction  layer  is  not  possible. 

In  Fig  13,  the  reaction  layer  thickness  is  calculated  and  illus¬ 
trated  based  on  the  results  given  by  the  model  for  the  reference 
parameters.  This  graph  clearly  shows  that  the  expression  given  in 
(72)  for  the  calculation  of  the  reaction  layer  thickness  is  relevant 
and  accurate. 

Two  important  comments  can  be  made  on  this  length.  Firstly,  it 
is  not  dependent  on  inlet  current  density  and  on  the  thickness  of 
the  electrode:  as  A  is  proportional  to  L2,  the  ratio  of  L  to  the  square 
root  of  A  is  independent  of  L,  which  confirms  the  results  given  in  Fig 
12.  Away  from  this  reaction  zone,  the  electrode  is  inactive  and  this 
gives  a  valuable  indication  for  the  choice  of  its  thickness.  Secondly, 
Lreac  depends  on  parameters  of  the  process  that  are  difficult  to 
obtain  experimentally  such  as  exchange  current  density,  i°,  or 
specific  reactive  volume,  ejp-  Therefore  experimental  determination 
of  this  length  can  help  to  determine  those  quantities.  Note  also  that 
a  quick  parametric  analysis  can  be  made  on  Lreao  as  A  directly  de¬ 
pends  on  the  parameters  of  the  process. 


x(10'4m) 

Fig.  12.  Comparison  between  the  two  models  on  transfer  current  density  along  the 
anodic  cermet  for  different  electrode  thickness.  X-axis:  position  in  the  anodic  cermet 
(on  the  left  the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  transfer  current 
density  multiplied  by  the  specific  reactive  volume. 


712 


M.  Dumortier  et  al.  /  Journal  of  Power  Sources  248  (2014)  703-713 


Reference  calculation 


A  =  324 
L=  150  jam 
L  =  58.3  jam 

reac  r 


-2.0 


-1.5 


Reaction  layer 


10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 


-1.0 


-0.5 


4° 


> 

3( 

u> 


x(10‘4m) 

Fig.  13.  Comparison  between  the  two  models  on  transfer  current  density  along  the 
anodic  cermet  for  different  conductivities  ratios.  X-axis:  position  in  the  anodic  cermet 
(on  the  left  the  anodic  channel,  on  the  right  the  electrolyte).  Y-axis:  transfer  current 
density  multiplied  by  the  specific  reactive  volume. 


4.  Comparison  with  results  from  the  literature 

The  1 D  model  was  compared  with  results  from  other  studies  on 
two  scales: 

•  The  macroscale  level  ( I—V  curves)  where  we  compared  the  I—V 
curves  given  by  the  1 D  model  with  experimental  and  correlated 
curves 

•  The  microscale  level,  where  we  compared  the  results  and 
models  in  literature  for  the  estimation  of  the  effective  electrode 
thickness. 


4.1.  Macroscale  level 


For  the  macroscale  level  comparison,  the  expression  of  potential 
losses  across  the  electrode  assembly  for  the  ID  model  is  required.  It 
is  obtained  by  adding  the  ohmic  losses  in  the  metallic  and  the 
ceramic  phase  by  integrating  Equations  (24)  and  (25)  using  Equa¬ 
tions  (58)  and  (59)  for  each  component  of  the  assembly  (anode, 
electrolyte,  cathode). 


/ 


ao 


(74) 


The  1 D  model  is  then  run  and  compared  with  the  experimental 
results,  using  the  data  given  by  Meng  Ni  et  al.  [30]  in  Fig  14  (for 
proton  conducting  SOECs)  and  Y.  Patcharavorachot  et  al.  [31]  (for 
proton-conducting  SOFC)  in  Fig  15.  The  experimental  conditions 
(i.e.  operating  conditions  and  electrodes  properties)  required  for 
the  calculation  are  usually  not  detailed  enough  on  experimental 
studies,  therefore  the  required  data  have  been  taken  from  the  both 
previous  studies  cited  above  [30,31]. 

The  ID  model  shows  good  comparison  with  the  experimental 
results  and  the  correlations.  On  Fig.  14,  the  experimental  results 
deviate  from  the  models  and  experimental  results  given  by  Meng 
Ni  et  al.  [30]  for  high  current.  However,  the  nonlinear  behaviour  of 
the  experimental  curve  could  be  caused  by  an  increased  hydration 
of  the  ceramic  due  to  higher  steam  fluxes  in  the  electrodes,  hence 
increasing  the  proton  conductivity  of  the  ceramic  and  therefore 
decreasing  the  ohmic  losses  in  the  ceramic.  It  is  not  the  case  in  the 
models,  where  this  proton  conductivity  is  kept  constant. 


Fig.  14.  Comparison  between  theoretical  simulation  results  and  experimental  data 
given  by  Ni  et  al.  [30]. 


4.2.  Microscale  level 

Several  experimental  results  and  models  can  be  found  in  liter¬ 
ature  on  the  effective  electrode  thickness  of  the  reaction  in  a  fuel/ 
electrolysis  cell  electrode.  The  definition  of  the  effective  thickness 
varies  from  a  study  to  another,  however  experimental  results 
usually  range  the  value  of  the  effective  thickness  between  a  few 
microns  to  a  few  hundreds  of  microns,  which  is  consistent  with  our 
study  (Lreac  =  58  pm).  Lu  et  al.  32]  and  Fukunaga  et  al.  [33]  reported 
better  cell  performance  with  larger  electrode  assuming  this  was 
caused  by  a  larger  effective  thickness.  In  our  model  though,  the 
effective  thickness  remains  independent  of  the  electrode  di¬ 
mensions.  It  is  not  obvious  to  know  which  observation  is  right:  this 
difference  may  either  results  from  the  deviation  of  the  1 D  model  at 
high  currents,  or  from  the  electrode  properties  varying  during  their 
preparation  prior  to  the  experiment.  In  addition,  the  correlation 
between  a  better  performance  and  a  larger  effective  thickness  is  not 


- Present  model  -  1073K 

a  Experiment  [32]  -  1073K 
—  Model  [32]-  1073K 

- Present  model  - 1 173K 

<  Experiment  [32]  -  1 173K 

. Model  [32]-  1173K 
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Fig.  15.  Comparison  between  theoretical  simulation  results  and  experimental  data 
given  by  Patcharavorachot  et  al.  [31]. 
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obvious.  Some  other  studies  also  observed  that  the  effective 
thickness  vary  with  the  value  of  the  inlet  current  density,  which 
again  may  be  the  case  for  high  currents,  where  the  ID  model  is  not 
valid  anymore. 

Models  in  literature  are  usually  based  on  equivalent  resistance 
modelling  and  present  different  expressions  of  the  effective  elec¬ 
trode  thickness.  These  models  are  usually  not  derived  from  a  set  of 
equations  and  are  built  to  extract  the  specific  resistances  of  the 
systems  by  fitting  experimental  data. 

Juhl  et  al.  presented  in  1998  a  resistance  model  used  for  estab¬ 
lishing  the  expression  of  the  effective  thickness  as  the  ratio  of  the 
triple  phase  boundary  resistance  to  the  resistance  of  the  ceramic 
[34].  The  physical  meaning  of  this  expression  is  similar  to  the 
physical  meaning  of  the  A  number  presented  in  this  study  in 
Equation  (62)  and  used  for  the  calculation  of  the  effective  thickness 
in  Equation  (73).  Juhl  et  al.  also  expected  the  effective  thickness  to 
increase  with  temperature,  which  is  the  case  in  our  model.  Shika- 
zono  et  al.  also  noted  that  the  effective  thickness  changed  with  gas 
composition,  which  is  consistent  with  our  study  as  well  35]. 
Finally,  O’Hayre  et  al.  [36]  and  Adler  et  al.  [37]  presented  models 
where  the  value  of  the  effective  thickness  decreases  with  an 
increasing  exchange  current  density,  which  again  is  consistent  with 
the  model  presented  in  this  work. 

It  is  however  difficult  to  assess  which  expression  is  the  most 
exact.  On  one  hand,  the  model  developed  here  was  rigorously 
derived  from  a  set  of  general  equations  for  low  currents  but  re¬ 
quires  a  lot  of  data  to  be  checked  experimentally.  On  the  other 
hand,  models  currently  reported  in  the  literature  are  probably  not 
comprehensive  enough  as  they  are  built  to  fit  the  results  given  by 
the  experiment  of  a  single  or  a  few  studies.  Although  all  of  these 
expressions  are  approximate,  including  the  one  presented  in  this 
study,  they  tend  to  show  the  same  dependencies,  which  is  a  good 
point.  In  addition  they  all  predict  values  of  effective  thicknesses  in 
the  same  range. 

5.  Conclusion 

The  analysis  presented  in  this  work  is  based  on  a  differential 
equation  system  developed  by  Dumortier  et  al.  [17]  to  model  the 
transport  processes  inside  solid  oxide  cells  using  a  cermet.  A  set  of 
assumptions  based  on  the  parameters  of  the  system  was  estab¬ 
lished  in  order  to  linearize  the  complex  differential  system  for  low 
inlet  current  densities.  By  integrating  the  system,  we  established 
explicit  functions  relating  the  quantities  of  the  process  to  the 
electrode  properties  and  the  position  in  the  cermet.  The  functions 
show  a  good  agreement  with  the  values  obtained  with  the  non- 
linearized  system.  Characteristic  dimensionless  numbers  were 
found  by  expressing  the  linearized  system  in  a  dimensionless  form. 
Together  with  the  system  equations,  we  described  a  relevant 
dimensionless  term,  A,  comparing  the  resistance  of  the  electrode  to 
charge  transfer  at  the  triple  phase  boundary  and  the  electrical 
resistance  of  the  metal  and  the  ceramic  to  charge  conduction.  This 
dimensionless  term  was  shown  to  be  pertinent  to  the  study  of  the 
distribution  of  the  electrochemical  reaction  rate  across  the  cermet. 
We  also  showed  that  the  electrochemical  reaction  can  be  confined 
to  a  layer  whose  dimension  is  calculated  with  this  number.  The 
thickness  of  this  layer  does  not  depend  on  inlet  current  density  or 
on  the  thickness  of  the  electrode.  It  is  therefore  an  indicator  to 
optimize  the  thickness  of  the  electrode.  In  addition,  the  layer  has 


been  explicitly  formulated  with  the  parameters  of  the  process  and 
the  properties  of  the  cermet.  Consequently,  the  values  of  quantities 
that  are  difficult  to  obtain  experimentally  can  be  determined  by 
experimental  measurement  of  the  layer.  A  comparison  with 
experimental  data  and  correlations  from  the  literature  on  this 
effective  thickness  showed  that  our  model  is  in  good  agreement 
with  the  current  results  on  the  matter. 
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